PARALLEL ALGORITHMS FOR SOME COMPUTATIONAL 

GEOMETRY PROBLEMS 


by 

A. Pavan Kumar Reddy 


<c t / TA 

ft a4<sP 

css 

i99Z 

& 

(j?E2 



DEPARTMENT OF COMPUTER SCIENCE Sc ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

FEBRUARY, 1995 



Parallel Algorithms for some 
Computational Geometry 
Problems 


A Thesis Sub'^i^e^ | ^ 
in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


by 

A. Pavan Kumar Reddy 


to the 

DEPARTMENT OF COMPUTER SCIENCE AND ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

Feb, 1995 



2 2 t'An |£95/c.5£ 

CENTRAL I '~9ARY 

I ! r KA 

<** a. .119115 


C_S£~ l°\<) r-fM- R££>- 



Acknowledgements 


I would like to express my sincere gratitude and 
heartfelt thanks to my thesis advisor. Dr. Sanjeev Saxena, for 
extending constant guidance through out this work. He introduced 
me to the interesting field of Parallel Algorithms. I am grateful 
to him for patiently explaining all my trivial doubts. I also 
learnt (at least tried) art of technical paper writing from him. 

I thank my examiners and Dr. Tapan K. Sengupta and Dr. Asish 
Mukhopadhyay for careful reading of thesis and for their 
constructive suggestions which has led to improved presentation 
of this thesis. I am also extremely grateful to Dr. Mukhopadhyay 
for 

a) Pointing out a special case (and how it should be handled) 
of one algorithm 

b) suggesting possibility of limitation (or error) in another 
algorithm 

c) drawing attention to a papers by Goodrich and by Chazelle 

I would like to thank Prof. S. Biswas for his excellent 
teaching. I am grateful to Dr. Mukhopadhyay and Dr. R.K. Ghosh 
from whom I learnt many things. 

I express my sincere thanks to all my friends at I.I.T., 
Kanpur for the delightful company they have given. I am thankful 
to Sudershan Banerjee for the useful discussions that I had with 
him. I thank Srin, RS, Parag and Singhai for helping me in need. 
As a prominent member of Gaali group I am thankful to my 
co-members TAN, TRK, Neti and Govind, for making me experience 
colours of life. I thank Balcha, Kathi, Kraja, Shankar and 
K(presi) for enriching my cultural life. 

I would like to express my deepest thanks to my parents, 
brother and all my relatives for the constant love that they have 
shown. I am very grateful to my cousin Pradeep, who has always 
been a constant source of inspiration for me. Finally I thank 
Lady Luck for favouring me. 


(A. Pavan Kumar Reddy) 



fi.ii npur. 

VW±& 


Certificate 


This is to certify that the work contained in the thesis titled Parallel algorithms 
for some computational geometry problems by Aduri Pavan Kumar Reddy 
(Roll No: 9311101), has been carried out under my supervision and that this work 
has not been submitted elsewhere for a degree. 


Feb 15, 1995 



5 — 


Dr. Sanjeev Saxena 
Assistant Professor, 
Department of Computer 
Science and Engineering 
UT, Kanpur 



1 


Abstract 

Parallel algorithms for following problems are proposed in this thesis. 

1. An optimal O(logn) time algorithm to find tangents from a set of points to a 
Convex hull is obtained, provided the point set satisfies some properties. 

2. It is shown that shortest path tree in a weak visible polygon can be computed 
in optimal O(logn) time with processors on a CREW PRAM without 
triangulating the polygon. 

3. An optimal parallel algorithm to triangulate a star shaped polygon is obtained 
on a CREW PRAM. It works in 0(log n) time with processors. 

4. A simple parallel algorithm to find convex hull of a simple polygon that runs 
in O(logn) time with processors is presented. Model of computation is 
CREW PRAM 

5. Some fast parallel algorithms are proposed to compute visibility polygon from 
a point. First a trivial 0(o;(n)) time algorithm with n 2 processors is described. 
Then a O(loglogn) time algorithm with processors is described. It is 
shown that this algorithms can be improved to get a algorithm that takes 
0((Iog Iogn) 2 ) time with processors where e is a constant (> 1). Later it is 
proved that this algorithm can be easily modified to obtain an optimal O(logn) 
time algorithm. All these algorithms run on COMMON CRCW PRAM 

6. The reverse problem of computing point set from a Voronoi diagram is consid- 
ered. A fully optimal parallel algorithm is obtained, on a CREW PRAM. 
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Chapter 1 


Introduction 


Computational geometry is concerned with the computational complexity of geomet- 
ric problems within the frame work of analysis of algorithms. Before the birth of 
Computational Geometry, mathematicians studied geometry with the intention of 
discovering properties of geometric objects. With the emergence of computers many 
researchers started working in the ares like pattern recognition, computer graphics 
and image processing where the basic stress was on the design of efficient algorithms 
based on the geometric properties of the objects involved. As classical treatment of 
geometric objects was not well suited to the design of efficient algorithms, it became 
necessary to select the useful concepts that are suitable for the purpose of the field 
of the algorithms. This marked the birth of discipline “Computational Geometry”. 

Many problems in this field need to be solved very fast, as these problems come 
from many real life applications. But we are now at the limits of what can be achieved 
through sequential computation. So emerged the field “Parallel Computational ge- 
ometry” to cope up with our needs. As many algorithms designed for sequential 
computer, are inherently sequential one needs to develop new paradigms suitable for 
parallel computation. 

In this thesis parallel algorithms for some visibility and shortest path problems, 
are studied. These problems have many applications in computer graphics, robotics 
and pattern recognition. Some of the applications are : 

• Suppose a robot is to be moved from one point s to another point t in a polygon. 
And the robot can not go out of the polygon. The robot should follow the 
shortest path between the two points with above constraint. This is equivalent 
to finding the shortest path in a polygon. 

• Suppose a communication channel is to be set up between two points in a simple 
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polygon whose sides are opaque to transmission. The transmission originating 
at a point s, called source is to be transmitted to a point t, called destination. 
Clearly, a direct communication is impossible unless the points s and t are in 
sight of each other inside the polygon. So some repeaters are to be installed 
in side the polygon to make communication feasible. What is the minimum 
number of repeaters required. This problem is equivalent to minimum link 
path between two points in a polygon. Minimum link path in a polygon is the 
path that has minimum number of vertices. 


1.1 Basic definitions 

In this thesis, the term “polygon” P, refers to the set of points lying on the boundary 
and inside it. Boundary of polygon P is denoted by bd(P). 

• Convex hull of a point set is the smallest convex polygon that includes all points 
of the set. 

• A polygon P is said to monotone with respect a line /, if any perpendicular line 
from / intersects bd(P) at most two places. 

• If the edges of a polygon do not intersect except at end points, and each vertex 
shares exactly two edges then it is called simple polygon . 

• Two points in a polygon are said to be visible from each other if and only if the 
line segment joining them completely lies inside the polygon. 

• A polygon is star shaped polygon if there exists a point c, inside it such that all 
points in P are visible from c (See Figure 1.1a). 

• Kernel of a polygon is set of all points from which the entire polygon is visible. 

• A point p is said to be weakly visible from an edge uv of polygon if there exists, 
a point c on uv such that p is visible from c. 

• The polygon P is called weak visible polygon if all points of P are weakly visible 
from an edge uv of polygon (See Figure 1.1b). 

• Shortest path in a polygon is the shortest Euclidian path which does not intersect 
any edge of polygon. 

• Triangulation of a polygon P, is division of P into maximum possible number 
of triangles, such that the vertices of triangles are vertices of P. If P has n 
vertices it can be divided into n — 2 triangles. 
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Figure 1.1: a: Star shaped polygon, b: Weak visible polygon 


• The expression right(pi,pj,p k ) ( left(pi,pj,p k )) is true if the point p k lies to 
the right (left) of the ray piPj. 

• Voronoi diagram of a point set S, is the partition of the plane into regions such 
that each region is the locus of points (x,y) closer to a point of S than any 
other point of S. 

1.2 Parallel computation models 

All sequential computers are based on the abstract model Random Access Machine 
proposed by von-Neumann. It has a read only input tape, on which the program is 
stored and write only tape on which the machine writes the output of the program, 
and memory. It stores intermediate values in memory. It can do operations like 
Addition, Subtraction, Comparison, Jump ( from one portion of program to another 
portion) on the contents of its memory. It can also read from input tape and write 
on output tape. The flow of the program in RAM model is sequential, i.e it can not 
execute more than one instruction at a time. 

A Parallel Random Access Machine employs p synchronous RAMs, all having unit 
time access to a shared memory. All the RAMs share the common program stored 
in read only tape. At any point of time the RAMs can either read from or write in 
the common memory or execute a step of the program. There are different types of 
PRAMs which differ in their capacity to do simultaneous reads and writes. From now 
we call a RAM as processor. 

An Exclusive Read Exclusive Write (EREW) PRAM allows neither simultaneous read 
nor simultaneous write by more than one processor to the same memory location. 

A Concurrent Read Exclusive Write (CREW) PRAM allows simultaneous reads but 



4 


not writes. 

A Concurrent Read Concurrent Write (CRCW) COMMON PRAM allows both si- 
multaneous reads and writes to same location by more than one processor, but a write 
attempt is valid only if all the processors try to write the same value in that location. 

If the worst-case running time of the best known sequential algorithm for a prob- 
lem is t, an optimal parallel algorithm for the same problem runs in £ time using p 
processors. An optimal fully parallel algorithm for the same problem runs in constant 
time with t processors. 


1.3 Overview of thesis 

In Chapter 2 the problem of finding tangents from a set of points to a convex hull is 
considered, when the point set has a special property. This problem seems to have 
many applications. We use this as sub-step in algorithms presented in Chapter 3 
and Chapter 4. An optimal O(logn) time algorithm with linear work is proposed 
for this problem. In Chapter 3 an optimal algorithm which uses processors 
and takes O(logre) time to computed shortest path tree in a weak visible polygon 
is described. An optimal O(logn) time algorithm for triangulating simple polygon is 
presented in Chapter 4. A simple parallel algorithm to find the convex hull of a simple 
polygon using the concepts of visibility is obtained in chapter 5. All algorithms in 
these chapters work on CREW PRAM. In chapter 6 some fast parallel algorithms to 
find the visibility polygon from a point are presented. First a trivial 0(a(n)) time 
algorithm is described. Then an O(loglogn) time algorithm with processors is 
discussed. It is also shown that visibility polygon from a point can be computed in 
O((loglog n) 2 ) time with processors. With a slight modification of this algorithm 
an optimal O(logn) time algorithm for the same problem is obtained. In chapter 7 a 
^ fully parallel algorithm for finding the point set from a Voronoi diagram is described. 



Chapter 2 


Finding Tangents from a Set of 
Points to a Convex Hull 


2.1 Introduction 

In this chapter we describe a problem which will be used as a sub-step in Chapter 3 
and Chapter 4. It is known that tangent to a convex hull from a point can be easily 
computed in O(logn) time. In fact tangent between two convex hulls can also be 
computed in (logn) time [29]. Here we consider the problem of finding tangents from 
all vertices of a polygon to a convex hull, where the vertices of the polygon have a 
special property. Assume that vertices of convex hull C are indexed ci,c 2 , ••*,c n in 
clockwise direction. Also assume that the vertices of P are indexed pi,p 2 , • • * ,p m 
in anti-clockwise direction. Chain(p,-,p,) denotes the boundary of P formed by the 
vertices Pi,P2, •* - ,Pm in anti-clockwise order. Chain(c,-, cj) is also defined similarly. 
Let the boundary of polygon be denoted by bd(P). In this chapter we first describe a 
sequential algorithms that runs in 0(m + n) time, then an optimal parallel algorithm 
which takes 0 ( log n) time is proposed. 


2.2 The algorithm 

The problem of finding tangents to a convex hull C from an arbitrary point set with 
linear cost is open. However, it is possible to obtain a linear time algorithm, if the 
vertices of the polygon satisfy following property. 
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Property 2.1 If tangent from a point p,- to C touches C at Cm then tangent from 
any point pj where j > i touches C at a point c n where n < m. 

First consider sequential version of the algorithm. 

1. Find the tangent from pi to C. Let it touches C at 

2.0. Assume that we have proceeded up to p,- and we have to find tangent 
from p,+j. Let the tangent from p,- touches C at Cj. 

2.1. If right (cj + i, Cj, p, + i),then Pi+\Cj is the tangent i=i+l, 

GoTo 2. 

2.2. If left(c J+ i, Cj, p,+i) then j=j+l, GoTo 2. 1 

It can be seen that the above algorithm works in 0(m + n) time 
Now consider the parallel version of the algorithm. 

1. Divide bd(P) into parts each having logn points. Call first point 
in each part tail and last point head. 

2. Find tangents to C from head and tail of each part. 

3. for all parts pardo 

3.1 Consider each part as Chain(p,-, pj) with p,- as tail and pj 

as head. Let tangents from p; and pj touch the convex 
hull C at c t - and cj respectively. Divide the chain chain(c,-, 
cj) into parts, where car(i,j ) is number of points 

in Chain(c,-, Cj). Define head and tail similarly for each 
part of C. 

3.2 For each head in Chain(c,-, Cj) find the point p* in Chain(p,-, 

Pj) such that the tangent from head passes through p*. 

REMARK: Thus, we have divided P into parts, each part having at most logn 
points. The tangents from heads of these parts divide C into parts but each part 
may have any number of points. Step 3.1 divides C into parts such that each 
part has at most logn points. Step 3.2 divides further divides each part of bd(P) 
such that each subpart has at most log n points, and overall P is divided into at most 
parts. C is also divided into at most parts. 

3.3 Take corresponding parts in P and C and apply sequential algorithm 
to each part. 

Claim 2.1 Above algorithm works in O(logn) time with processors. 
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Proof: Step 2 can be done in O(logra) time with m/log n processors. Step 3. 2 can 
be done in O(logn) time with processors. In this step for a point c,- € C we have 
to find a point from a group of at most logn points, which can be done in O(logn) 
time with a single processor. So over all algorithm works with above processor and 
time bounds. I 



Chapter 3 

Shortest Path Tree in a Weak 
Visible Polygon 


3.1 Introduction 

Finding shortest paths in a polygon is one of the basic problems in Computational 
Geometry. Shortest path problems arise in many application areas such as computer 
graphics, robotics and vision. Routine for shortest path is also used in other compu- 
tational problems like, finding the visibility information [16], and minimum link path 
in a polygon [34]. Therefore, a great deal of research work has been devoted to finding 
efficient algorithms for shortest path and its variants [16, 34] (like link paths and L\ 
paths). 

Lee and Preparata [26] describe a linear time algorithm to find shortest path 
between two vertices in a simple polygon. Guibas et a/. [23] has shown that the shortest 
path tree rooted at any point in a polygon can be computed in linear time. Both 
algorithms assume that polygon is triangulated. As Chazelle [11] recently proposed a 
linear time algorithm for triangulating a simple polygon we can sequentially compute 
shortest path tree rooted at any vertex in linear time. Ghosh et al. [18] describe a 
linear time algorithm for finding shortest path tree in a weak visible polygon, without 
triangulating the polygon. They solve the problem with a new characterisation of 
weak visible polygon in terms of shortest path between two arbitrary vertices of the 
polygon. 

Attempts have also been made to obtain parallel algorithms for these problems. 
Goodrich et al. [22] obtain a parallel algorithm for computing the shortest path tree in 
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O(logn) time with processors(on CREW model), if the polygon is triangulated. 
But the best known parallel algorithm for triangulating a polygon takes O(logn) 
time with n processors [21]. As an optimal algorithm for triangulating weak visibility 
polygon is not known, we cannot use the algorithm of Goodrich et al. to get an 
optimal algorithm even if the polygon is weak visible. Here we propose an optimal 
O(logn) time algorithm for finding shortest path tree in a weak visible polygon. Our 
model of computation is Concurrent Read Exclusive Write (CREW). We use the 
characterisation of Ghosh et al. to solve our problem. The algorithm uses rootish 
divide and conquer paradigm [5, 20] to reduce resource bounds. 

We assume that the weak visible polygon P is given as (array of) counterclockwise 
sequence of vertices (i.e., their respective x and y coordinates) p 1} p 2 , • • 'i Pn • Assume 
that the polygon is weakly visible from the edge pip 2 . The symbol P denotes both 
the boundary of the polygon (bd( P)) as well as the region bounded by it. Chain(p,-, 
Pj) denotes the boundary formed by the vertices p,-, p, + x, • • •, pj in counter clockwise 
order. The shortest path between two vertices of polygon should not intersect any 
edge of the polygon. The shortest path between two vertices p t - and pj is denoted as 
SP(pi,pj). Given any three points p< = (x,-,t/;),pj = (xj,yj) and p k = (x k ,y k ), let 
S — Xk(yi — i ij) + y k (xj — Xi) + yjXi — yiij. If S' < 0 then we say p k lies to right of the 
line passing through p,- and pj. If S > 0, then it lies to the left of the same line. If 
the interior angle at a vertex p,- is convex then u; is called a convex vertex. An edge 
PiPi+i is called convex edge if both p,- and p,- + i are convex vertices. Two vertices are 
visible from each if and only if the segment joining them lies completely inside P. A 
point p is said to be weakly visible from an edge st, if there is a point z on the line 
segment st such that p and z are visible. If every point in P is weakly visible from st 
then P is said to be weakly visible from st and P is called a weak visible polygon. 


3.2 Properties of weak visible polygons 

We assume that the edge p\p-i is convex edge. If a vertex p,- of P is weakly visible 
from convex edge pjp 2 of P, then the following properties hold [23] 

Property 3.1 1.1 5P(p 2 ,pJ makes right turn at every vertex in the path 

1.2 SP(pi,pi) makes left turn at every vertex in the path 

1.3 SP(pi,pi) and 5P(p 2 ,p,j are two disjoint paths and they meet only atpi 
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1.4 The region enclosed by SP(pi,Pi) and SP(p 2 ,Pi) and p\p<i is totally contained 
inside P. 

Ghosh et al. [18] prove following additional properties of weak visible polygon. 
Property 3.2 2.1 For any two vertices pi andpj, no vertex , except those of Chain (pi^pj) 

can be on SP(pi,Pj). 

2.2 For any two vertices pi and pj, SP(pi,pj) makes right turn at every vertex in 
the path. 

Thus from property 3.2, for any two vertices pi and pj, SP (pi,Pj) is the same as 
the convex hull of Chain(p,', pj). 

Let T{ be the shortest path tree of polygon chain- Chain(p r ,p r+m ), rooted at 
vertex p T (for any vertex p r ). Further, let pi be any internal node in this tree, with 
children pj,p?,* •• ,p{- We define an ordering on these children as follows: draw a 
straight line from p, to any child p*. If x is the number of children of p; lying to 
right of the line p.pf then we define rank of pf to be x. We also define an ordering 
among leaf nodes. Rank of a leaf node U is zero if every node in the path from p r to 
li has a rank zero. Let /,• and lj be any two leaf nodes, with lowest common ancestor, 

L = LCA(U,lj). Let us assume that the rank of first node in the path from L to lj 
is greater than the rank of first node in the path from L to f;. Then rank of /,• and 
lj are defined in such a manner that rank of lj is greater than rank of /;; note that 
ranks are in range 0 to number of leaves in tree —1. 

If pi is a leaf node with rank / then nextleaf(pi) will denote the leaf node with 
rank / + 1. If p t - is an internal node then SuccT(pi) will denote the child of p,- 
which has rank zero. For a leaf node p;, SuccT(pi) is the first node in the path from 
lowest common ancestor of: p; and nextleaf(pi), to nextleaf(pi). And 5ucc(p,) = 
the successor of pi in P (— p,+i) for all nodes. In 2] the leaf with rank zero is called 
the left most leaf and the leaf with maximum rank is called the right most leaf. 

Lemma 3.1 For any node p,- in tree T{, SuccT(pi) = Succ(pi). 

Proof: First consider the case when p; is not a leaf node. We prove the Lemma by 
induction on the level of node p,-. 

Base case: The node p,- is at level zero, i.e it is root p r . By Property 3.2, shortest 
path from p,- to p,+i (= Succ(pi)) should be the edge p,-p;+i. Thus Succ(pi) is at level 
one. Assume that Succ(pi) is not 5uccT(p,), but some other vertex pj. Note that 
there is an edge p r -ip r in the polygon. Now there should be a polygon chain from pj 
to p r _ i (to form the polygon P), and the tree edge p,SuccT(p,) should lie inside P. 
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Figure 3.1: Proof of Lemma 3.1 


Further it should be possible to reach p r _ x from pj along this chain during anti-clock 
wise traversal. Observe that the chain from pj to p r _ 1 can be formed in two different 
ways (Figure 3.1). In the first case, the tree edge piSuccT(pi) lies outside P\ and 
PiSuccT(pi) can not be an edge of shortest path tree, a contradiction. 

In the second case the tree edge p{SuccT(pi) lies inside P, but in this case we are 
doing a clockwise traversal to reach p r _i, which means that the polygon vertices are 
indexed in clockwise order, this contradicts our basic assumption that P is given as 
anti-clockwise sequence of vertices. Hence SuccT(pi) = 5ucc(p ,‘ ) when p,- is root node. 
This completes the base case. 

Induction argument: Let the claim be true for all nodes up to level h. Let p,- 
be any node at level at h + 1. Assume that Succ(pi) is not SuccTfa), but some 
other vertex pj. This implies that the vertices in path from p r to p,- and then to. 
Pj form a polygon chain. So shortest path from p r to pj should pass through the 
vertices of this chain only. But it is passing through some other vertices (vertices in 
the tree path from p r to pj); again a contradiction. So pj cannot be Succ(p,), and 
Succ(pi) = SuccT(pi). 

The proof when p; is a leaf node is similar to that of base case. I 

Corollary 3.1 In a tree T the vertex with maximum index is the right most leaf. I 
Note that if we know the shortest path from p; to pj and also the shortest path from 
Pj to Pkt where both are convex paths, then to find the shortest path from p< to p* 
we have to merge the two hulls by finding the tangent. 
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Figure 3.2: The sub-problem, Finding tangents to C,- 


Finding shortest path tree of the given polygon requires solution of following 
intermediate problem: 

Problem 3.1 Given a convex hull <7,-, and a layer of hulls (here layer of hulls is the 
tree Ti) find tangents from each point in the layer of hulls to <7,- .(see Figure 3.2). 

We solve the Problem 3.1, using processors in logn time, where m is total 
number of points, in section 3.3. 

Lemma 3.2 Assume that the polygon P is divided into k polygon chains, of con- 
secutive vertices. Suppose that shortest path tree Ti of each of chain i is known 
(i = 1, • • • , k), and that the right most leaf of tree T;_i is the same as the root of tree 
Ti. Assume that root of first tree Ti is p r . Then shortest path tree of the polygon P 
rooted at p r can be obtained by merging trees T\, 7*2, *• • , T* in O(logn) time, with 
max(P, processors on a CREW PRAM; here n is the total number of vertices 
in polygon P. 

Proof: Let r,- be the root of TJ. Observe that we can look each tree T; as layers 

of hulls (see Figure 3.3). First we will merge the upper most hulls of all the trees 
(i.e. first we will find the shortest path from p T to some vertices that lie on the upper 
hull of each tree.); for this we may consider only the upper most hull and ignore all 
other vertices in the tree. Assign one processor to each pair ( TfTj ) for j < i and find 
tangents from Ti to all Tjs, where j < i. Consider the set of tangents to T;. Divide 
these tangents into two sets. First set consists of tangents that pass through some 
vertex p,- (of T<). Second set consists of tangents that pass through the vertices pj 
where j < i. Discard all the tangents of second set. Of all the tangents in first set 
keep the tangent which makes maximum angle with x-axis. Discard other tangents 
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Figure 3.3: Merging shortest path trees 

to T{. This gives a partial shortest path tree. Now we have to find shortest path from 
P T to the vertices in the lower hulls (of T{) and to the vertices that lie between r t - and 
Pi on the upper most hull, of T; (for all trees). Consider (7,-, the shortest path from 
p r to r;_x where r t _i is root of T< (or right most leaf of 7)_j). To find shortest path 
p r to all other vertices of T t , it is sufficient to find tangents from all these vertices of 
Ti to C{. As we can solve Problem 3.1 in same resource bounds, we are done. I 
REMARK: By Corollary 3.1, the right most leaf of tree T;_i is the same as the root of 
tree Ti. 

Thus, we can find the combined shortest path tree of the fc-chains in 0 ( log n) time 
with max(fc 2 ,m/logn) processors, where m is total number points in all the trees. 
The algorithm to compute shortest path tree in P is: 

1. If the problem size is less than logn then solve sequentially, else, divide the 

polygon chain into parts, each part having at most ^/n\ogn points. 

2. Solve each sub-problem recursively 

3. Merge all yj%^ sub-problems, as described in Lemma 3.2 
Claim 3.1 Above algorithm works in C(logn) time with processors. 

Proof: By Lemma 3.2, the sub-problems can be merged in O(logn) time. The 

recurrence equations are 

T(log n) = log n 

T(n) = T(\/n logn) -f logn 

Solving the recurrence the equation we get T(n) = O(logn). To merge k sub- 
problems we need max(& 2 ,m/ logn) processors (Lemma 3.2). As k is always less 
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Figure 3.4: P 1 P 2 is not convex edge 

than yjn log n and m is less than n, processors are sufficient. I 

So far we have assumed that edge pip? is convex. Let us relax this condition. If 
P1P2 is not convex, then extend the edge pip? so that it intersects bd(p) at p a and 
Pb as shown in Figure 3.4. Delete chains Chain(p 2 , Pa) and Chain(p{,,pi) and find 
the Shortest path tree in the new polygon. We have to find shortest paths to the 
vertices in the chains, Chain(p 2 ,p 0 ) and Chain(p(,,pi). Note that the polygon formed 
by Chain(p 2 , Pa) and the edge p 0 p 2 forms a star shaped polygon with p 2 as a star 
point. Similarly the polygon formed by the chain Chain(pi ( ,p 1 ) and the edge pxPb 
forms a star shaped polygon with pj as star point. We can easily triangulate these 
star shaped polygons as the star point is one of the vertices and find shortest path 
tree in these polygons [22]; basically, we will be done if we join vertex p 2 with all 
other vertices in the polygon. 

Although, the algorithm as described, computes the shortest path tree rooted at 
p 2 , we can easily modify the algorithm to compute the shortest path tree rooted at 
any vertex. Let us assume that we have to compute the shortest path tree rooted 
at vertex p t . Consider all vertices that lie in the chain, Chain(p,-,p 1 ) and find the 
shortest path tree rooted at p,- for these vertices as above. To find the shortest paths 
for the vertices that lie in the chain Chain(p 2 ,p,), reverse the order of vertices so that 
now they are arranged in clockwise order. Use the same algorithm, except that roles 
of “left” and “right” are interchanged. Thus we have the following Theorem: 
Theorem 3.1 In a weak visible polygon the shortest path tree rooted at any vertex 
can be found in O(logn) time with processors. 
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Figure 3.5: Proof of Property 3.3 

3.3 Solution of Problem 3.1 

In this section we consider the solution of Problem 3.1. We are given a convex hull 
Ci, and a layer of hulls (here layer of hulls is the tree T{) and are required to find 
tangents from each point in the layer of hulls to Ci (see Figure 3.2). 

Let us assume that tree T,- has m hulls (i.e T, has m leaves). We will divide the 
points in all layers (of hulls) into two sets. If p t is top most point of a hull ( of Ti) 
then we put the points from r,- to p t in the first set, and rest of the points of the hull 
in the second set. We describe the algorithm to find tangents from points in the first 
set; tangents from points in the second set can be similarly found. We do not know 
of any linear time sequential algorithm for finding tangents from an arbitrary set of 
points to a convex hull. However the points in the problem under consideration are 
vertices of a weak visible polygon, using Property 3.3, it is possible to find tangents 
with linear cost. 

Property 3.3 If tangent from a point p,- to Ci touches Q at p m then tangent from 
any point pj where j > i should touch Ci at a point p n where n <m. 

Proof: Let the tangent from p,- to Ci touch Ci at p m . Let pj ( j > i ) be any other 
point. Note that tangent from p, can pass through point pk (k > m) (see Figure 
3.5) only if pj lies to left of tangent p;p m . In this case it can be easily observed that 
there exists an edge p c p<* ( c, d > i) which intersects the tangent p,p TO (Figure 3.5). 
This implies that shortest path from p r to p, is passing through some other vertex 
which does not lie in the polygon chain Chain(p r ,p,), which contradicts Property 3.2 
of weak visible polygon. I 

Note that this property is the same as Property 2,1 of Chapter 2. Hence Prob- 
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lem 3.1 can be solved optimally. This completes the proof of Theorem 3.1. 



Chapter 4 

Triangulation of Star Shaped 
Polygons 


4.1 Introduction 

Algorithm for triangulating simple polygons finds applications in many other prob- 
lems in Computation Geometry. The algorithm for triangulation is used as a pre- 
processing step in algorithms like planar point location [24, 27, 29]. Triangulation 
algorithm is also used in finding visibility information and computing shortest paths 
[16, 26]. Clearly, these problems can be solved more efficiently, if the input is given 
as a triangulated polygon. On sequential computational model, many algorithms are 
available to triangulate a simple polygon. Garey et al. obtain an O(nlogn) time 
algorithm for this problem. They solve the problem by dividing the polygon into 
monotone polygons, and then triangulate them. Fournier and Montuno [15] also de- 
scribe an algorithm with same time complexity; here, the polygon is divided into 
trapezoids and then into one sided monotone polygons. Tarjan and Van Wyk [35f 
first break the 0(n log n) bound by obtaining an 0(n log log n) time algorithm to tri- 
angulate a polygon. They also obtain a linear time algorithm if the polygon is weakly 
visible from an edge. Recently Chazelle [11] has obtained a linear time algorithm to 
triangulate a simple polygon, thus solving one of the well known open problems in 
computational geometry. If the polygon has special structure then it possible to ob- 
tain simpler algorithms (with smaller constant factors); extremely simple linear time 
algorithms exist for triangulation, when the polygon is monotone or star shaped [15]. 

In parallel setting, Goodrich [21] obtains an O(logn) time algorithm with n pro- 
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cessors to triangulate a polygon, on a CREW PRAM. The algorithm works even when 
polygon has holes. If the trapezoidal decomposition of polygon is given then the al- 
gorithm takes, O(logn) time with j— processors. Yap [38] has also obtained another 
algorithm with same resource bounds. If the polygon is monotone then the problem 
can be solved optimally in O(loglogn) time [8] on a CRCW PRAM. In this chapter 
we propose an optimal algorithm for triangulation, if the polygon is star shaped. Our 
model of computation is the CREW PRAM. 

We assume that the star shaped polygon P is given as (array of) counterclockwise 
sequence of vertices (i.e., their respective x and y coordinates) pi, p 2 , •• p»; z(p«) 

and y(pi) denote x and y co-ordinates of the point p;. Assume symbol P denotes both 
the boundary of the polygon ( bd(P )) as well as the region bounded by it. Chain(p,-, 
Pj) denotes the boundary formed by the vertices p,-, p t+ i , • • -, pj in counterclockwise 
order. CH(p,-,pj) stands for the convex hull of Chain(p,-,pj). Two vertices are visible 
from each other if and only if the segment joining them lies completely inside P. The 
polygon P is called star shaped if there exist a point c inside P such that the entire 
polygon P is visible from c; c is called kernel of the polygon. If a horizontal line is 
drawn through a vertex p,-, and if adjacent edges of p,- lie on both sides of the horizontal 
line then p,- is called a regular vertex . The vertex p,- is called stalagmitic vertex if 
both adjacent edges lie below the horizontal line, and vertex p,- is called stalactitic 
vertex if both adjacent edges lie above the horizontal line. We call a vertex p,- bottom 
peak vertex if it is a stalactitic vertex and it lies below the horizontal line through 
the kernel point c (these definitions are from [15]). We call top peak vertex if it lies 
above the kernel point c and it is a stalagmitic vertex. 


4.2 The algorithm 

We triangulate a star shaped polygon by decomposing the polygon into smaller and 
more simpler objects like monotone polygons. 

Without loss of generality, assume that the left most point of the polygon is p\. 
Draw a horizontal through the kernel point c. Let it intersects the edge p a -iPa on 
the left side and the edge pt-iPb on the right side (note that the line should intersect 
the polygon at only two points as P is star shaped). We call Chain(p 0 ,pj) the lower 
chain and the Chain(p(,,p 0 ) the upper chain (see Figure 4.1). 

Fact 4.1 If pj and pj are any two vertices of l.ower chain such that t < j then the 
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Figure 4.1: The star shaped polygon P, and lower hull 

convex hull of the vertices that lie in Chain (p;,Pj), (CH(p,-, pj)^does not intersect any 
other part of the polygon P. 

Proof: The kernel point c can at most touch the CH(p,-,p ; ) but cannot lie inside 
it as Chain(pi>?i) * s P ar t °f the lower chain. Now assume that there is some edge 
(PiP ) the polygon P that intersects an edge of CH(p l -,p J ). This means that one 
of the vertices pi or p m is not visible from c, a contradiction. I 

Consider the bottom peaks of the lower chain, let pit, be the bottom peak with 
minimum index. Let pu be the top peak with maximum index in upper chain. Then 
observe that 

Fact 4.2 The chain, Chain(pn,p«>) is monotone with respect to y axis. 

Proof: All vertices in the chain, Chain(p/i,pjfc) are regular vertices. This is so 

because if a vertex is not regular then it should either be a stalagmitic vertex or a 
stalactitic vertex; this implies that the vertex will be either a top peak, or a bottom 

peak. 

By definition of regular vertex, if p, is a regular vertex then y(p,_i) < j/(p,-) < 



Similarly we can prove that Chain(p r i,p rt ) is also monotone with respect to y 
axis where p r b i s k he vertex with maximum index among bottom peaks and p rt is the 
vertex with minimum index among top peaks. 

We next describe the algorithm. 

Step 1 Find upper and lower chains, by finding vertices p„ and pt,. 

REMARK: This can be done in O(logn) time. 
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Step 2 Find the hull of lower chain (portion lying between p a and pt,). The hull 
divides polygon P into two parts the upper part (P u ) and the lower part (Pj). 
We will consider the problem of triangulating P/; P u can be triangulated in a 
similar manner. 

Observe that if the portion of hull between p a and pt has k edges, then polygon 
Pj is in turn divided into k sub polygons (Pj , • •• ,P/ fc ). Each sub polygon is 
bounded by the hull edge and the polygon chain lying below that edge. If the 
polygon Pj has X; points then assign processors to it. 

REMARK: As in this step we are finding the convex hull of a simple polygon it can be 

implemented in 0(log n) time with processors [37] 

Step 3 For all polygons pj, • • • ,pf do 

Step 3(a) Identify all non-regular vertices in Pj and find non-regular 
vertices having the maximum ( p r ) and minimum index ( pi ). Also 
find the vertex pj having the minimum index among all vertices of 
Pj and the vertex p e having the maximum index among all points 
of Pj. 

Step 3(b) Find following convex hulls — CH(p/,pj), the convex hull 
of Chain(p/,p/); CH(p;,p r ) the convex hull of Chain(p/,p r ) and 
CH(p r ,p c ) the convex hull of Chain(p r ,p 4 ). 

Step 3(c) Triangulate following polygons — 

the polygon bounded by chains Chain(p/,p r ) and CH(p/,p r ); 
details of triangulation algorithm will be described 
later 

the monotone polygon bounded by chains Chain(p/,p/) 
and CH(p/,p/) (see Fact 4.3 below) and 
the monotone polygon bounded by chains Chain(p,.,p 4 ) 
and CH(p r ,p e ) (see Fact 4.3 below) 

Monotone polygons can be easily triangulated [21], 

Step 3(d) In O(logn) time, find tangents from p t , the top most point 
of CH(pj,p r ), to both CH(p/,p/) and CH(p r ,p e ). 

REMARK: We are left with the polygon as shown Figure 4.2b. 


Step 4 Triangulate remaining polygon. 
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(c) 

Figure 4.2: Steps during execution of algorithm, (a) The polygon P[ after step 3(c); 
(6) After step 3(d); (c) Monoconvex polygon 

REMARK: We show how to triangulate remaining polygon later. 

Fact 4.3 In polygon P/, the chains, Chain(p/,p/) and Chain(p r ,p e ) are monotone 
with respect to y axis. 

Proof: The proof is similar to proof of Fact 4.2. I 

We next consider implementation of Step 3(b), the problem of triangulating P,tai> 
the polygon bounded by Chain(p/,p r ) and CH(p/,p r ). Actually P ata i is divided into 
smaller polygons by hull edges of CH(p/,p r ). Call the part of P tta i that lies below any 
hull edge as P JU (>»ta/ (See Figure 4.3a for an example of P, u b,t a i)- Next we look at the 
problem of triangulating P lu bstai (implementation of Step 3(b)) 

3(b). 1 Identify all bottom peaks {pbpu" • ,Pbpr) of P»ub.t«\- Find Pa,Pt2," * ,Ptr, top 

most points (i.e stalagmitic point) between two consecutive bottom peaks. If 

y(pti- 1 ) > y{Pti) then draw a horizontal line from p t i so that it touches the 

Chain(pti_i,p6 pt ) at (say) p mp ,-. Similarly if y(p t i) < y(pti+i), then draw a line 

to touch the Chain(p 6p , +1 , p ti+1 ) at p mpi+ i C.NTK n L LI BRAR> 

l t. f., KANPU R 

Z a HUB 
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3(b). 2 Triangulate each polygon that lies between the horizontal segment p m ipu and 
Chain(p mt -,pt,). 

3(b) .3 Join two successive topstalagmitic vertices; a vertex p t i is a topstalagmitic ver- 
tex, if y(p ti -i) < y(p ti ) > y(pu+i) 

3(b). 4 Triangulate each P S3U bstah the sub-polygon that lies below two successive top- 
stalagmitic vertices. The sub-polygon is a monotone polygon, (See Fact 4.4, 
below). 

3(b). 5 Triangulate remaining monotone polygon 

It can be easily shown that polygon in Step 3(b).2, the polygon lying between 
horizontal segment p m ipu and Chain(p m ,-,p tl ) is monotone and hence, can be easily 
triangulated. The remaining polygon after Step 3(b). 2 is as shown in Figure 4.3b. 
Let pu and p t j (*' < j) be two successive topstalagmitic vertices. To triangulate each 
Pssubstah (to implement Step 3(b).4), observe that 

Fact 4.4 The chain formed by vertices pu , p m ,+i , p ( ,+i , • • • , Ptj is combination of at 
most two monotone (with respect y axis) chains. 

Proof: We show that in this chain (say TSChain(p t ,-,p tJ )) the y co-ordinate first 
decreases and then increases. If p t k is any vertex in this chain then observe that 
the vertices between p t k and p m k+i are regular. Consider any three successive top- 
stalagmitic vertices p t k-i,Ptk and p t k+i, then either y(p t k-i) < y(ptk) < y(ptk+ i) or 
y(Ptk-i) > y(ptk ) > y(ptk+i)- 

Note that if y(p t k- i) < y(Ptk) and y(p t k) > y(Ptk+i), then p t k becomes a topsta- 
lagmitic vertex, and this vertex marks the end of chain TSChain(p tI -,p t j). Thus, the 
chain TSChain(p t ;,p tJ -) can be a combination of at most two monotone chains. I 

Thus, the polygon bounded by the line pti,ptj and TSChain(p t ,-,p tJ ) is a monotone 
polygon, and can be easily triangulated; this completes the description of Step 3(b).4 
The remaining polygon (after Step 3(b).4) is monotone as the chain formed by joining 
two successive topstalagmitics is monotone. As Step 3(b). 5 can be done in optimal 
O(logn) time, entire Step 3(b) can be implemented in O(logn) time with pro- 
cessors. 

Coming to implementation of Step 4, let the tangent from the top point (say 
p t ) of CH(p/,p r ) touch hulls CH(p/,p/) and CH(p r ,p e ) at p x and p v respectively. 
As sub-chain Chain(p/,p r ) of chain CH(p/,pj) is monotone, and also as sub-chain 
Chain(p„,p e ) of chain CH(p r ,p c ) (by Fact 4.3) are monotone, and the polygon formed 
by Chain(p/,p I ) the line p E p { , the line p t p y and the Chain(p v ,p e ) is also monotone 
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and can be easily triangulated. 

We have to triangulate the polygon that is bounded by Chain(p x ,p/) of CH(p/,p/)and 
Chain(p;,p t ) of CH(p/,p r ) and the tangent p t p x . The polygon has two convex chains 
(call them the left and right chains) sharing a common vertex, and a straight line 
joining the other ends of the chains. We call such polygons monoconvex polygons 
(see Figure 4.2). Similarly the polygon bounded by Chain(p t ,p r ) of CH(p t ,p r )and 
Chain(p r ,p y ) of CH(p r ,p y ) and the tangent p t p z is a monoconvex polygon. Monocon- 
vex polygons can be triangulated as follows: 

4.1 Find tangents from each point of right chain to the left chain. 

4.2 Let tangents from the points p,- and p; + i pass through pj and p* respectively. 

Then Join the point p t - with all the vertices that lie between pj and p*. 

Observe that the vertices of right chain satisfy Property 2.1 (Chapter 2), as the right 
chain is a convex chain. This implies that Step 4.1 can be optimally implemented 
using the algorithm for finding tangents from a set of points to a convex hull (see 
Chapter 2). 

Hence the following theorem can be stated: 

Theorem 4.1 A star shaped polygon can be triangulated in O(logn) time with 
processors on a CREW PRAM. 




Chapter 5 


Convex Hull of a Simple Polygon 

5.1 Introduction 

Finding convex hull of a set of points is a well studied problem in Computational 
Geometry. Several algorithms are available to solve this problem efficiently on a 
sequential computational model [29], As sorting is linear time reducible to convex 
hull problem convex hull problem has a lower bound of fl(nlogn). Moreover if the 
points are vertices of a simple polygon then the problem becomes much simpler, and 
can be solved in linear time [29]. Convex hull can also be computed in linear time, 
in case points are sorted (with respect to one of the co-ordinates) or are vertices of a 
monotone chain. 

On parallel computational models also, the problem of finding convex hull has 
been extensively studied [2]. Chow gave an algorithm that works in O(log 2 n) time 
with n processors. Atallah and Goodrich [5] have obtained an optimal 0(log n) time 
algorithm on CREW model, to find the convex hull of a set of points in the plane. 
They use rootish divide and conquer technique to bring down the time complexity.- 
Aggarwal et al. [1] also describe an optimal algorithm for the same problem. Cole 
and Goodrich [13] propose an optimal algorithm for the same problem, without using 
the square root function. They use cascading divide and conquer paradigm [4] for 
merging of sub-problems in O(logn) time. An optimal O(logn) time algorithm with 
processors is given by Goodrich [20] to compute the hull when the points are 
sorted. He brings down the processors by stopping the recursion after sometime and 
then solving the sub problems sequentially. The algorithm also works when the points 
are vertices of a monotone chain. Miller and Stout [28] obtain an optimal O(logn) 


25 



26 


time algorithm on EREW model, to compute the convex hull. The problem of finding 
convex hull of a simple polygon is considered by Wagner [19, 37] and he proposes an 
optimal O(logn) time algorithm for it. Here in this chapter we obtain the same 
bounds, but the technique used is different. We use the concepts of visibility to solve 
the problem. Atallah et al. [3] describe an optimal O(logn) time algorithm with 
processors on EREW PRAM model, for finding visibility polygon from a point. 
Another result needed is an optimal O(logn) time algorithm with linear cost (on 
CREW) by Goodrich et al. [22] for finding shortest path in a triangulated polygon. 
In this chapter we use CREW PRAM model of computation. 

We assume that the (simple) polygon P is given as counter clockwise sequence 
of vertices (i.e., their respective x and y coordinates) pi, p 2 , •••, p„. The symbol P 
denotes both the boundary of the polygon (bd(P)) as well as the region bounded by 
it. Two points in P are visible from each other if and only if the segment joining 
them lies completely inside P. The visibility polygon of P from a point c is that part 
of p which is visible from c. Let cp; be a tangent from a point c to convex hull of 
P. The segment cp, is called right(left) tangent of P if all points of P lie to the right 
(left) of the line passing through c and p,-. Unless stated otherwise in the chapter, line 
passing through p, and pj direction will be from p, to pj. Chain(p,-, pj) denotes the 
boundary formed by the vertices p,-, p, +1 , •• •, pj in counter clockwise order. UH(P) 
will denote the upper hull of polygon P. The left most and right most points of P 
will be denoted by pi and p r respectively. 


5.2 The algorithm 

If C is convex hull ( with vertices v\, • • • ,v m in counter clockwise direction) and v is 
a point outside it then observe that: 

Fact 5.1 If the left and right tangents to C from t; passes through V[ and v T then all 
the vertices in Chain(u/,v r ) of C are visible from v. 

Fact 5.2 If the left (right)tangent from v to C passes through vi (tv) then v lies to 
the right(left) of vi-\vi (v r -itv) and to the left (right)of (tvv r +i) 

Now we describe the algorithm to compute the convex hull of P. 

1. Find pi and p r , the left most and right most points of P. 

2. Consider all the points p { of P. Compute the anti-clockwise angle (clockwise) 
angle between lines pipi and pip r . Let p c (p«j)be the point for which the angle is 
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(c) 

Figure 5.1: (a): The polygon P; (b): The polygon P'; (c): The polygon p". 
maximum. 

3. Let c be a point that lies to left of ptp c and to the right of p r Pd- Join c with 
Pi and p r forming a new polygon (P') c,p/,p/ +1 , • • • ,p r ,c. Find the visibility 
polygon P"of P' from c 

4. Triangulate the polygon P" . 

5. Find the shortest path from p a to p\> in P", which is the upper hull. 

First observe that 

Claim 5.1 If a point p, belongs to the upper hull of P then it is visible from c. 
Proof: The points pi and p T are present in the upper hull of P. Let C be the convex 
polygon formed by UH(P) and the line segment p T pi. The point c lies to left of the 
line passing through p\ and p c , it also lies to the right the line passing through p T pi- 
Thus due to Fact 5.2 tangent to C from c passes through pi. Similarly we can see that 
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the right tangent from c to C passes through p r . Consider the polygon P c formed by 
cpt, UH(P), and p T c. No other point of P can properly lie ( i.e. it should not even 
touch the edges) inside P c , because if it lies then it should belong to UH(P). Now 
there are no objects (polygon edges) in side P c . Tangents from c pass to C through 
pi and p r . By Fact 5.1 all points in the chain Chain(p/,p r ) of C are visible from c. 
Thus all the points in UH(P) are visible from c. I 

Theorem 5.1 Convex hull of P can be computed in 0(log n) time with proces- 
sors. 

Proof: First we prove the correctness of the algorithm. By claim5.1 upper hull 

of P is same as upper hull of P" (formed during step 3). So we will show that for 
the polygon P" (formed in step 3) convex hull and shortest path from pi to p r axe 
same. Observe that if pipj is an edge of the hull UH(P) then pi is visible to Pj in the 
polygon P'. So if p,pj is an edge of UH(P) then shortest path from p,- to pj in P" is 
nothing but the edge PiPj. We prove that if p,pj and pjPk are two consecutive edges 
of UH(P) then shortest path from p,- to Pk in P" is p,- — Pj — p*. First we prove that 
Pk is not visible from p,- in P". Draw a straight line from p,- to p*, note that p* should 
lie above this line, as p,- ,pj and p* are consecutive edges of convex hull. Now there 
is a chain, Chain(p,-, pj) (of P") below the hull edge p.pj, similarly there is a chain 
Chain(pj,pfc) (of P") below the hull edge PjPk- So the horizontal line segment from 
Pi to Pk intersects both these chains, so p,- is not visible from pk- This implies that 
the shortest path from p; to pk can not lies below p,pj nor can it lie below PjPk- So 
it should lie above (or on) these edges. This means the shortest path from p,- to pjt 
should lie in P". But in P" there are no vertices of P other than the vertices that 
are in UH(P). So the shortest path from p,- to Pk should pass through the vertices of 
UH(P) only. As UH(P) is convex chain the shortest path from p; to p*. should pass, 
through pj. So the shortest path between p,- and pk is pi - pj - pk, which is equal to 
convex hull of Chain(Pj, Pj) (of P"). With similar arguments we can show that the 
shortest between pi and p r is equal to convex hull of Chain(p/,p r ) of P". This proves 
the correctness of the algorithm. We shall now see the implementation details. 

Steps 1 and 2 can be done trivially in O(logn) time with processors. Step 3 
can be done using the algorithm of [3] with ^ processors in O(logn) time. In step 4 
we have to triangulate the polygon P". This is a star shaped polygon as we have 
found the visibility polygon of P' from c to form P" and the star point is c. So we 
can triangulate it easily, we have to just join all the vertices of P" with c. We can 
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find the shortest path in a triangulated polygon in O(logn) time optimally [22] So 
the entire algorithm works in 0(log n) time with processors. I 



Chapter 6 

Visibility Polygon from a Point 


6.1 Introduction 

Finding visible regions of a polygon from a source is well studied problem in com- 
putational geometry. Several linear algorithms have been obtained for finding visible 
regions of a simple polygon from a point source or from a line source. Lee [25] shows 
that visibility polygon from a point can be computed in linear sequential time. Avis 
and Toussaint [6] described a linear time algorithm to test whether a given polygon 
is weakly visible from a given edge of the polygon. Suri and Sack [31] obtain a linear 
time algorithm to test whether a given polygon is weakly visible from any edge of the 
polygon. Guibas et al. [23] proposed linear time algorithm for finding weak visibility 
polygon from a segment inside polygon. Ghosh [16] extends the notion of visibility 
from a point or line segment to visibility from convex set, to get linear time algorithms 
for computing complete, strong and weak visible polygons from a convex set. 

Attempts have also been made to obtain parallel algorithms for these problems. 
Atallah et al. [3] obtain an optimal 0(log n) time algorithm for finding visibility poly- 
gon from a point on the Exclusive Read Exclusive Write (EREW) model. Chandru 
et al. [10] propose an O(log(m + n)) time algorithm with (m 4- n) processors on the 
CREW model for finding complete visibility polygon. 

Here we give some fast parallel algorithms for computing visibility polygon of a 
simple polygon from a point. We first describe an 0(a(n)) time algorithm with n 2 
processors. Then we propose a fast algorithm for computing the visibility polygon in 
O(loglog n) time. This algorithm uses ^ processors. We propose another algorithm 
which takes 0((log log n) 2 ) time with processors. We also describe an optimal 
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O(logn) time algorithm. 

We assume that the (simple) polygon P is given as (array of) counter clockwise 
sequence of vertices (i.e., their respective x and y coordinates) pi,p 2 , • • • ,p n . The 
symbol P denotes both the boundary of the polygon (bd(P)) as well as the region 
bounded by it. Two points in P are visible from each if and only if the segment joining 
them lies completely inside P. If c is a point in P then the visibility polygon of P from c 
is the set of all points in P that are visible from c. Chain(p,-, pj) denotes the boundary 
formed by the vertices p,-,p,+i, • • • ,Pj in anti-clockwise order. UChain(p,-,p,) denotes 
the boundary formed by the vertices between p; and pj] the direction can be either 
clockwise or anti-clockwise. In this chapter, the point c will be assumed to lie inside 
the polygon P. 

In Section 6.2 we describe the 0(a(n)) time algorithm. In the Section 6.3 we 
propose an O(loglogn) time algorithm that uses processors. In Section 6.4 we 
describe an algorithm which uses processors and takes a time of O((loglogn) 2 ). 
In the same section we show that this algorithm can be easily modified to give an 
optimal 0(log n) time algorithm. 

6.2 0(a(n)) time algorithm 

First we consider the following “trivial” algorithm. The algorithm takes 0(a(n)) time 
with n 2 processors on a CRCW PRAM. 

Let c be the point from which we are computing the visibility polygon. 

0. Initialise the visibility polygon P' to be P. 

1. For each vertex p,- of the polygon P, check whether the segment cp,- 
intersects any edge of the polygon, if cp,- intersects an edge then mark 
Pi as ‘invisible’ else ‘visible’. 

2. For each vertex p, find the (first) vertex pj (j > i) such that pj is visible 
from c, and no other vertex in Chain(p,-,pj) is visible from c. 

3. Find the intersection of the edge Pj-iPj with the segment cp,-. Let pf-_ x 
be the intersection point; in P\ replace Pj_i with p' t and mark j/j_ x 
‘visible’. 

4. The polygon formed by joining adjacent ‘visible’ vertices of P' gives 
the visibility polygon of P from c. 

Now we will show that the algorithm correctly computes the visibility polygon 
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from c. In first step we have identified the vertices of P that are visible from c. 
Consider any vertex p, that is marked ‘visible’, let pj be the first ‘visible’ vertex that 
we encounter in counter clockwise walk from p,-. Now we have to compute the visible 
part of Chain(p,-,pj). We know that no polygon vertex that is in Chain(p,-,pj) is 
visible from c. Because if any vertex in that chain is visible from c then p, cannot 
be the first visible vertex encountered during the anti-clockwise walk from p,-. So 
the visible part of Chain(p,-,pj) is just the visible part of edge Pj-iPj- Hence, it is 
sufficient to compute only the intersection of edge Pj-iPj with cpj. Thus, the polygon 
formed by visible vertices gives the visibility polygon of P from c. 

To analyse the algorithm, assume that we have n 2 processors at our disposal. 
Assign n processors to each vertex of the polygon. We can test whether a point pi is 
visible from c in constant time using n processors assigned to p,-; basically, we have to 
test whether the segment cp,- intersects any side of the polygon. As there are n edges 
in polygon, we can test in constant time. Step 2 can clearly be done in constant time. 
All vertices marked ‘visible’ in the output array, give the desired polygon. To delete 
vertices that are not visible from c, for each visible vertex we have to find the nearest 
left and nearest right vertices that are marked ‘visible’. This is equivalent to the 
problem of finding nearest zeros in a binary array and can be solved in 0(a(n)) time 
either using the recursive star data structure of Berkman and Vishkin [9] or using 
the algorithm of Ragde [30]. Hence, we can find the visibility polygon of P from any 
point c in 0(a(n)) time using n 2 processors. 

6.3 The O(loglogra) time algorithm 

Assume that we have divided the polygon into z sub-chains Pi, • • • , P z . Let P{ be the. 
visibility chain (part of chain visible from c) of P t - from c; in definition of P[, we do 
not consider other chains Pj, j ^ i (i.e., each P, is not looked upon as an obstacle, 
but is assumed to be transparent). Let pi be the first and p r be the last point (end 
point) of the chain, P- respectively. 

Let P" be the polygon formed by adding segments (c,p/) and (c,p r ) to P(. Consider 
the intersection of polygon P/' with chains Pj’s (1 < j < z, and j ^ »). Consider 
the parts of these chains that lies inside P/' (they can touch the edges of P[')\ we call 
these parts Subp,- chains (here Pi indicates that these chains are inside P/’, similarly 
we can have Subpjchains). We will divide these Subp, chains into three types based 
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(«) ( 6 ) 

Figure 6.1: a : Types of Subp, -chains, b : Forming P"' 


on their intersection type with P" (see Figure 6.1a). 

Type 1: The Subp, chain touches segment cp/. 

Type 2: The Subp^chain lies inside polygon P". It touches neither cp r nor cp/. 

Type 3: The Subpichain touches cp r . 

Let pj ( pi ) be any vertex that belongs to a Type 3 (Type 1) chain; and the line 
cpj (cp,)makes an angle (anti-clockwise direction) of Angj (Angi) with cp/. Find the 
vertex p m i n (Pm&x) such that the value of Ang„^ n (Ang^x) is minimum (maximum) 
among all Ang^s {AngC s). 

Let the line cp„m X (cp m in) intersects P{ at pj e /t (Pri 3 ht)- Note that the vertices 
in Chain(p/p max ) and also the vertices in Chain(pnu n p r ) are not visible from c (Fig- 
ure 6.1b). Let the polygon formed by cp^, Chain(p m ,^ p m j„) and p m ;„ c be P"'. 

Consider any Subp/chain that lies inside polygon P"'. Let x (y) be the point 
of Subpichain such that cx (cy) makes minimum (maximum) anti-clockwise angle 
with cpt among all other vertices belonging to the same chain. Call this minimum 
(maximum) angle as StartAngle (EndAngle). If the StartAngle of one Subpichain is 
equal to the EndAngle of another Subpchain then treat them as a single chain, and 
call it ContinuousSubpichain. 

Figure 6.2a gives an example of ContinuousSubpichain. In the figure there sure two 
Subpichains Chain(p,q) and Chain(r,s). Here StartAngle of one chain is equal to En- 
dAngle of another chain. So we should treat both of them as a single ContinuousSubpichain. 

Lemma 0.1 There can be at most two ContinuousSubp.chains in P-". 
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Figure 6.2: a : Example of ContinuousSubpidiains. b : Proof of Lemma 6.1 

Proof: We will prove the Lemma by contradiction. Assume that there are three 
ContinuousSubp.-chains (call them ChA, ChB and ChC) in P/" (Figure 6.2b). There 
should be a polygon chain from one end point of ChA to ChB, but it cannot com- 
pletely lie inside P/"; as in that case C hA and ChB will be in the same ContinuousSubp;chain. 
So let us assume that the chain from ChA to ChB does not completely lie in P{" 
(Figure 6.2b), In this case we can merge ChB, ChC and Chain(c,<f) to form a 
ContinuousSubp, chain. Hence, in both cases we cannot have three Continuous S ubp,chai ns . 

A similar proof can be given for non existence of more than three ContinuousSubp, -chains. 
Hence, the Lemma follows. I 

By Lemma 6.1, to find the visible portion of Pi (or P/) in the polygon P it is 
sufficient to find the end points of the two ContinuousSubpjchains (in case there jure 
two chains, otherwise one). Thus, we can find the region not visible from c (these 
chains block some part of P,) and hence the region of P, visible from c. Less informally, 
the algorithm is: 

1. Divide the polygon chain into parts each having at most logn points and 
assign n processors to each chain P,. 

2. For each chain P,- find the visibility chain P{ of P, from c. While computing 
visibility chain of P; ignore all other chains (i.e. consider edges of other chains 
as transparent). 

REMARK: As each chain P, has at most logn points, P(, the visibility chain of 
Pi can be computed in 0(a(n)) time using the algorithm, of previous section. 

3. For each chain P/, join both its end points to c to form polygon P”. Find the 
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intersection of ff with all other chains. 

REMARK: Now we have Subp, chains for each P( 

4. Find polygon P[" (as defined in proof of Lemma 6.1) for each i. 

REMARK: Subp, chains can be either Type 1, Type 2 or Type 3. In case of 
Type 1 and Type 3, Subp, chains, we can form polygon P-" by finding vertices 
(say x and y) that make minimum (among Type 3 vertices) and maximum 
(among Type 1 vertices) angles with line cp/. Find intersection of cx and cy 
with P/ to get vertices p^x and p m j n . We can find the intersection of these lines 
with P( (in constant time) as there axe at most log n points in P- and we have 
n processors. This gives the desired polygon P"'. Type 2 Subp, chains will be 
inside the polygon P"\ 

5. Identify the ContinuousSubp, chains and find visible regions of Pi in P. 

6. Club the remaining regions of P( to form the visibility polygon of P from c. 

Claim 6.1 Above algorithm can be implemented in O(loglogn) time with pro- 
cessors. 

Proof: Step 1 can clearly be done in constant time with n processors. In Step 3, 
we have to find intersection of P" with all other chains. There are j— chains, and 
each P'i has n processors, thus we assign logn processors to each pair (p'/,p'-). 

We find intersection between each pair with log n processors as follows. Recall, 
that all chains are star shaped (as each chain is visible from c), and each chain has at 
most logn points. So the problems reduces to that of finding the intersection between 
two star shaped polygons each having at most logn points with logn processors. 
Maheswari and Ghosh [17] have shown that the problem of finding intersection of two 
star shaped polygons each having at most m points can be solved, in O(log m) time 
with processors. Thus, the problem of finding the intersection between each pair 
takes O(log log n) time. As we have to find intersection between P" with all other 
chains, all intersections can be found in O(loglogn) time with n processors. Thus, 
for all Pf s Step 3 can be done with ^ processors in O(loglogn) time. 

In Step 4, the functions minimum and maximum can be computed in O (log logn) 
time on CRCW model. 

In Step 5 we have to find the visible region of P( in P (considering the existence of 
all other chains; that is we assume that they are “opaque”). We will do this without 
explicitly computing the ContinuousSubp^chains. In P'" there can be at most 
Subp, chains. We have logn processors assigned to each Subp, chain. Let Chi be any 
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Subpichain with end points r\ and s,- (cr^ makes smaller angle than csi with the line 
cpi). Find the intersections tzj and y,- of cr^ and c$,- with P[. From now on we treat 
tij and Vi also as vertices of P/. Mark all the vertices (including ii; and u,-) and edges 
of P '■ as “visible”. Assume that with each vertex we can store two more data items, 
left item and right item. Now make Chi write right in “slot” for right item of u,- and 
left in the “slot” for left item of y,-. All concurrent writes are of “COMMON” type. 
Let ChA( with end points ChAf and ChA e ) and ChB (with end points ChBj and 
ChB e ) be two ContinuousSubpichains. 

Assume that cChA e (cCkBj) makes larger (smaller) angle than cChAj ( cChB e ) 
with the line cp/. Let the intersection of cChA e ( cChBj ) with P{ be Pa,- (P6,). Note 
that the right item of Pa,' will be empty, while the left item is filled with left. Also 
note that by Lemma 6.1 no other vertex can have this property. Similarly left item 
of Pbi will be empty, while right item is filled with right. And again by Lemma 6.1, 
this also cannot happen for any other vertex. So the visible portion of in P is the 
chain between Pa,- and p6,-. 

To implement Step 6, we delete the invisible parts of chain P/. The successor of 
Pa< will be ChA e and the predecessor of P6,- will be ChBj. To delete invisible parts 
efficiently we maintain chain P- as a binary search tree. Chain P( can be maintained 
as a binary search tree because vertices of P- are angularly sorted about c. Again by 
Lemma 6.1, we have to delete at most two sub-chains of P/. These trees are exactly 
identical to the hull trees maintained by Goodrich [20] and deletions and merging of 
trees can be carried as in [20]. 1 

Thus, the following theorem follows: 

Theorem 6.1 The visibility polygon of P from a point c lying inside P can be found 
in O(loglogn) time, with processors on CRCW PRAM. 

6.4 Reducing the number of processors. 

In this section we describe two algorithms, the first algorithm takes 0((log log n) 2 ) 
time with processors. Second is an optimal 0(log n) time algorithm. 

1 + 1 

6.4.1 The algorithm with processors 

In this section we reduce number of processors to (e is a constant > 1) with 
slight increase in time to 0((log log n) 2 ). Observe that in the previous algorithm 
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the costliest step is finding the intersection of two star shaped polygons. Here we 
will exploit a new property of these star shaped polygons to find their intersection 
efficiently; we use the fact that these star shaped chains are sub-chains of the original 
polygon P. 

Let us draw a horizontal line through c. Let p] (and p' T ) be the point of P/ such 
that the line cp\ (respectively line cp' r ) makes the minimum (maximum) anticlock wise 
angle with the horizontal through c, among all other points of p\. Similarly, we can 
define pj and pj for chain Pj. 

Lemma 6.2 Let P, and Pj be two chains of polygon P. Let P( and Pj be visibility 
chains of P, and Pj respectively. Then chain Pj can intersect chain P/ at most twice. 

Proof: We prove the Lemma by contradiction. Assume that Pj intersects P/ at 
three points (Figure 6.3a). Let the edges of Pj that intersect P/ be p a Pb, p c Pd and 
p e p/. As P is a simple polygon, none of these edge can be the edges of P. These 
edges were added when we computed the visibility chain of Pj from c. This implies 
that there exist a polygon chain (say) UChain(p a ,p(>) from p a to p&; similarly there 
exist two other chains UChain(p c ,pj) and UChain(p e ,p/). 

In first case, the UChain(p c ,p<j) is as shown (Figure 6.3a). In this case, polygon 
chain Pj should start at p c ; otherwise the chain, Chain(p s ,p c ) will be a part of Pj and 
will obstructs UChain (p/,pj); as a result, the edge p e pj will not be there. So we can 
assume that Pj starts at p c . 

Traverse polygon P in anti-clockwise direction, starting from p c . We will first 
encounter p* (and not p/, p<*, p*, or p a ). Now there should exist a chain joining p e and 
p/ (UChain(p e ,p/)), and the chain should intersect either cp t or cp r (Figure 6.3a); 
observe that this chain should be a sub-chain of Pj. In either case Pj cannot be visi- 
bility chain of Pj from c; because if UChain(p e ,p/) intersects cpi then UChain(pj,p/) 
is not visible; and if it intersects cp r then the UChain(pj,p(,) is not visible. So either 
edge p e pf or edge p a pb of Pj cannot exist. 

Similarly, in the other case, when the UChain(p c ,p (i ) is as shown in Figure 6.3b, 
we can prove that the edge p a Pb cannot exist. So Pj cannot intersect P/ three times. 
Similar arguments hold when the number of intersections are more than three. So Pj 
and P[ can intersect at most at two points. I 

In subsequent discussion we assume with out loss of generality that index p\. is 
greater then pj (otherwise interchange the role of “minimum index” and “maximum 
index” in following discussion). Let the angles made by the lines cp}, cp}., cpj and cpj 
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with the horizontal line be a, /?, 7, and 6 . The intersection between Pj and Pj can 
be any of the following six cases based on the relations between these angles. 

Case 1 (7 < a < /? < 0): There are three sub-cases as the chains may 
or may not intersect (see Figure 6.4c). 

Case 2 (7 < ot < 9 < /?): There are two sub-cases, as the chains may 
or may not intersect (see Figure 6.4a). 

Case 3 (7 < 6 < or < /?): In this Case the chains do not intersect (see 
Figure 6.4e) 

Case 4 (a < 7 < (3 < 9 ): Here there can be two sub-cases, as the 
chains may or may not intersect. This is symmetric to Case 2 (see 
Figure 6.4b). 

Case 5 (0 < 7 < $ < /?): Here also there are two sub-cases, as the 
chains may or may not intersect (see Figure 6.4d) 

Case 6 (a < < 7 < 9 ): Here the chains do not intersect (see 

Figure 6.4e). 

Observe that there cannot be any other case (or sub-case). For each case we show 
how to find the intersection between Pj and Pj. 

Cases 3,6 : Here there is no intersection. 

Case 5 : Here we have to find the intersection point of the line cpj with the chain 
Pj. If the intersection does not lie on the segment cp\ then Pj and Pj do not 
intersect, otherwise Pj lies in Pj'. 

Case 1 : In this case Pj is not visible from c or the intersection is nil; The case 
when the intersection is nil can be easily detected with binary search. If the 
intersection is not nil then we prove that P,- is not visible from c as follows: 

If Pj is as in Figure 6.6a then Pi is obviously not visible from c. However, if Pj 
is as in Figure 6.6b, then Pj intersects Pj at two points. Let the edges of Pj that 
intersect pj be p,„iPm2 and Pi n zPinA\ there should be a path from p' r to Pj. Let the 
path touch Pj at point (say) p^ tart . Let pJ tort lie in between vertices pj and p,„ 1 of 
Pj. Now the original polygon chain of Pj should leave chain UChain(pj,p,„i) of 
Pj at some point and it should enter either UChain(p, n 2,p >n3 ) or UChain(p in4 ,p^) 
(both chains are sub-chains of Pj). Let us assume that the original polygon chain 
of Pj first enter UChain(p; n2 ,Pm3) of Pj (see Figure 6.6b), and then leave the 
chain and then enter UChain(p, n4 ,pf) of Pj. Finally, the chain should reach p} 
to complete the polygon. The chain will reach p] as shown in Figure 6.6b; thus 



39 


it can be easily observed that P( is not visible from c. 

Note that if chain is as shown in Figure 6.6c, then c will lie outside P. 
Similarly it can be shown that P- is not visible from c in all other remaining 
cases. 

Cases 2,4 : As Case 2 and Case 4 are symmetric, we consider only Case 4. First 
find the intersection point of line cp\ with Pj, if the intersection point doesnot 
lie on the segment cp\ then the two chains do not intersect. Intersection of line 
cp' r with Pj, can be found by binary search; if we have processors then 

binary search takes O(loglogn) time (car(Pj) is number points of Pj). If line 
cp' T does not intersect Pj, then P- and Pj do not intersect. Otherwise, we have 
to find the edge PiniPini of Pj that intersects P/. It is sufficient to identify either 
point p,„j or point p,„ 2 . This is done using Claim 6.2 (described below). 

Claim 6.2 In Pj, the visibility chain of Pj, if neither p, nl nor p,- n2 have the minimum 
index amongst the vertices of Pj then, the intersection between P/ and Pj is not 
required in computation of visible part of P/. 

Proof: Let us assume that neither p tnl nor p, n2 have the minimum index. There 
should be a polygon path from p’ to the vertex with minimum index in P- (say Pmin)- 
Note that Pmin cannot lie between p* and p, n2 ; as in that case Pj will not be the 
visibility chain of Pj (see Figure 6.5b). In case pmin lies in the polygon chain between 
pj and pi n2 (see Figure 6.5a), it can be easily seen that there exists another chain CH 
(say) such that the region blocked by chain CH includes the region (i.e., is a super 
set of the region) blocked by Pj. Thus, we need not compute the intersection of P( 
and Pj. ■ 

This completes the description of Case 4. Case 2 is similar. 

Algorithm to find visibility polygon of P from c is: 

1. Divide the polygon chain into n« parts, each part having at most n 1- « points. 
Assign processors to each part P,- 

2. If the number points in P, is greater than log log n then recursively find P/, the 
visibility chain of P;. Otherwise solve the problem sequentially. 

3. To merge visibility chains find the visible portion of each P{ in P 

as follows. As there are processors assigned to each P/, we distribute these 
processors as follows: processors are assigned to each pair (P/, Pj). Form 

the polygon P- 1 and find the intersection of P- 1 with all the Pj’s as described 
above. 
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4. Find the visible region of PlinP 

5. Merge all remaining chains 

Claim 0.3 Above algorithm can be implemented in O((log logn) 2 ) time with 
processors. 

Proof: Step 1 can trivially be done in required resource bounds. Intersection 

between P" and P*, in Step 3 is found by finding angles a, /3, 7 and 6, and identifying 
the intersection types as explained above. As already seen the intersection can be 
found using binary search in O(Ioglogn) time with ™*( car (^)’ car ( f V)> processors. 

Each P( and Pj has at most n i- « points and there are processors assigned 
to the pair (P-,Pj). Thus, Step 3 can be implemented in O(loglogn) time. 

Step 4 is identical to Step 5 of algorithm described in Section 3, and Step 5 of 
this algorithm is again identical to Step 6 of the algorithm described in Section 3. As 
already seen, both steps can be implemented in O(loglogn) time. 

Thus the recurrence equation is 
T(n) = P(n 1- <) -f log log n 
T(loglogn) = log log n 

Hence, 2» = 0((log log n) 2 ). ■ 

Putting everything together, we have the following theorem: 

Theorem 6.2 Visibility polygon of P from c can be computed in O((loglog n) 2 ) time 
with processors. 

6.4.2 Optimal O(logn) time algorithm 

We next modify 0((loglogn) 2 ) time algorithm to get an optimal O(logn) time algo- 
rithm with processors. 

1. Divide the polygon chain into parts , each part having at most ^/(n logn) 

points. Assign processors to each P. 

2. If the number points in P,- is greater than logn then recursively find P-, the 
visibility chain of Pi . Otherwise solve the problem sequentially. 

3. To merge visibility chains P[, P 2 ', • • • find the visible portion of each 

P( in P as follows. As there are yj[ processors assigned to each P/, we 
distribute these processors as follows: one processor is assigned to each pair 
(P/, Pj). Form the polygon P-' and find the intersection of P" with all the Pj’s 
as described above. 

4. Find the visible region of P{ in P 
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5. Merge all remaining chains 

Claim 6.4 Above algorithm can be implemented in O(logn) time with proces- 
sors. 

Proof: All the steps in this algorithm are identical to those in the O((log log n) 2 ) 
time algorithm, except here Step 3 (binary search) which takes <3 (log y/n log n) time 
as we have (only) one processor for each pair. 

It can be easily observed that remaining steps can also be executed in 0(log y/n log n) 
time. So the recurrence equation is 
T(n) = T(y/n\ogn) + log(V'nlogn) 

T(logn) = logrt 

Or, T(n) — O(logn). ■ 

Thus, the following theorem 

Theorem 6.3 Visibility polygon of P from c can be computed in 0 ( log n ) time with 
processors. 








Chapter 7 

Getting Co-ordinates of Points 
from a Voronoi Diagram 


7.1 Introduction 

Suppose that we are given a Voronoi Diagram G of a point set S. In this chapter, 
we consider the “reverse” problem of getting the co-ordinates of the point set S from 
the Voronoi Diagram. This problem may arise when the point set S is lost, but the 
Voronoi Diagram remains, and we have to regenerate the point set. 

Voronoi diagram of a point set S, is the partition of the plane into regions such 
that each region is the locus of points (x,y) closer to a point of S than any other 
point of S. Each point of the set S is associated with a region (say), Voronoi region 
of the diagram. To be more formal Voronoi region associated with point p, (say 
V(pi)) is intersections of half planes H(puPj ) (1 < j < n and j ± t), where H(pi,pj) 
is the half plane containing p,- that is defined by the perpendicular bisector of line 
PiPj. Voronoi diagrams are used to solve proximity problems like nearest- neighbor 
search [29]. Given a point set S its Voronoi diagram can be sequentially computed in - 
0(n log n) time [29]. The best known parallel algorithm takes 0(logn log logn) time 
with n ^giogn processors, and <9 (log 2 n) time with processors [14]. The vertices 
of the Voronoi diagram are called Voronoi vertices and the points of the point set are 
referred as Voronoi sites. A Region of the Voronoi diagram is also called face. 

In some cases it may not be possible to store both the point set and its Voronoi 
diagram due to memory limitations. This means we have to store only the point set 
and compute the Voronoi diagram whenever needed. But it is not very attractive 
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solution as computing Voronoi diagram is quite expensive. In this chapter we provide 
alternate solution to this problem; store only the Voronoi diagram and compute the 
point set when needed. We propose a fully parallel algorithm to compute the point set 
from the Voronoi diagram. The algorithm runs in 0(1) time with n processors on a 
CRLW PRAM. This model allows simultaneous reads but does not allow simultaneous 
writes. As the Voronoi diagram G is a planar graph, we assume that the planar 
embedding of Voronoi diagram is given, i.e. for each Voronoi vertex vj its adjacent 
vertices are given in counter-clockwise order. We are also given the co-ordinates of 
the Voronoi vertices. This definition of planar embedding is consistent with [12]. 


7.2 Computing the point set 

The following properties of Voronoi diagram of a point set S are known [29] (it is 
assumed that points are in general position — no 4 points are co-circular). 

1. Every Voronoi vertex is intersection of exactly three edges of the diagram, i.e. 
Voronoi diagram is regular of degree three 

2. Every Voronoi vertex (u) is center of the circle(say C(v) defined by three points 
of S 

3. For a Voronoi vertex v, the circle C(v) does not contain any other points of s 
Observe that if the Voronoi diagram G is having only two faces, i.e. S has only two 

points, then the Voronoi diagram G is just a straight line l (Voronoi edge), and the 
point set S is obviously not unique (any pair of points for which / is a perpendicular 
bisector will do). 

Consider the case when G has three faces, then G will be as shown below (see 
Figure 7.1); G has only one Voronoi vertex (v) and three Voronoi edges incident at. 
v. Let us put the origin of the co-ordinate system at v. If S = {p,q,h}, then p, q 
and h should lie on a circle centred at v; let the radius of this circle be r. Further, let 
the coordinate of p, q and h, be (r cos o, r sin a), (r cos 0, r sin 0) and (r cos 7, r sin 7) 
respectively. Let mi, m% and m 3 be the slopes of Voronoi edges, i.e., the equations 
of line corresponding to Voronoi Edges are y = m x x, y = m 2 x, and y = m 3 x. (See 
Figure 7.1). As each Voronoi Edge is perpendicular to the line joining two points of 

S, • A , 

sin a — sinp _ —1 
cos a - cos 0 
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<«) 


(b) 


Figure 7 . 1 : The Voronoi diagram 


sin fi — sin 7 
cos /? - cos 7 


and 


sin 7 — sin a 
cos 7 — cos a 

As each Voronoi edge is equidistant from (i.e. 
points of 5, 


-1 
m 2 

-1 

m 3 

is a perpendicular bisector of) two 


sin a + sin /? = mi (cos a + cos /?) 


and 


sin ft + sin 7 = m x (cos (3 + cos 7) 


sin 7 + sin a = m x (cos 7 -f cos a) 


These six equations can be solved to get a, ft and 7. However, the radius r is not 
determined by these equations. Thus, when the number of Voronoi regions is three, 
then also the set S is not unique (see Figure 7 . 1 ), as we can choose any r, but a, (3 
and 7 are fixed. Our main result is the following theorem: 

Theorem 7.1 A Voronoi diagram G uniquely determines the point set S such that 
G is a Voronoi Diagram of S, provided that G has at least four regions, or equivalently 
S has at least four points. 
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Further co-ordinates of points in the set S can be obtained in 0(1) time with n 
processors on a Concurrent Read Exclusive Write (CREW) model. 

Proof: Let us first consider the case when there are four Voronoi regions in G (see 
Figure T.l); here there are two Voronoi vertices, and there is a Voronoi edge joining 
them. Let the two Voronoi vertices be a and b. Let p and q be the points in the regions 
adjacent to edge ah. Now with respect to Voronoi vertex a, we can determine the 
angles a, 0 and 7; hence points p and q are determined except for radius r. Similarly, 
we can also determine the angles of points p and q — a', 0' and 7' with respect to 
Voronoi point b, also; again the radius r' has not been fixed by above equations. As 
p (and similarly q) should be the intersection point ray coming out of a and b (as 
angles are fixed, so are these rays), point p is unique and can be easily determined; 
similarly point q is unique and can be determined. Thus, the all points in set S can 
be uniquely determined from the Voronoi diagram G. In case, there are more than 4 
regions, we can repeat the procedure for each pair of adjacent regions, and determine 
co-ordinates of all points. 

As Voronoi diagram is a planar graph, number of edges e is of the same order as n, 
the number of Voronoi points. Thus, using 0(n) processors, we can then assign one 
processor to each Voronoi edge a 6, and solve the equations in 0(1) time to determine 
the co-ordinates of 4 site points (each Voronoi point is defined by 3 points, here two 
points p and q are common). It is quite possible that co-ordinates of some point are 
determined more than once. At each Voronoi vertex we store the co-ordinates of the 
respective sites. So at each Voronoi point at most 3 Voronoi sites are stored. But 
each site is stored at many Voronoi vertices; a Voronoi site p is stored at all Voronoi 
vertices which belong to V(p) (the Voronoi region of p). We show how to copy these 
Voronoi sites into an output array of size n, without concurrent writes. We use Euler 
tour technique [36]. For each Voronoi region (V(p)) we can find the edges that bound 
V(p), using slight modification of Euler tour technique [7, 33], as follows. 

for each edge (i j) pardo next.onJace[(ij)]:= next[reverse[(i,j)]] ( =next[(j,i)] ); 

As input is given such that for each Voronoi vertex we know the Voronoi vertices 
adjacent to it, this algorithm can be applied. This gives a linked list (say L(p)) for 
each region (V'(p)), such that the vertices of the linked list L(p) are incident on the 
region V(p). Now the head of the list L(p) writes the co-ordinates of p in its slot. 
This avoids concurrent writes as only one processor writes in each location. This 
completes the proof of the theorem. I 



Chapter 8 

Concluding Remarks 


The problem of finding tangents to a convex hull from a set of points considered in 
Chapter 2, seems to have some interesting applications. In this thesis the point set 
satisfies a special property. It needs to be explored whether this property can be 
relaxed. In particular, if we can find tangents from vertices of a star-shaped polygon 
to a convex hull inside the polygon in optimal 0(log n ) time, then it may be possible 
to obtain an optimal algorithm to find complete visibility polygon from a convex 
set [16]. 

In Chapter 3 an optimal parallel algorithm is described to compute the shortest 
path tree in a weak visible polygon. This algorithm does not use parallel triangulation 
algorithm and is optimal. It is still not known whether shortest path tree in simple 
polygon can be found without triangulating the polygon. Even all the sequential 
algorithms proposed first triangulate the polygon. To compute shortest path tree 
in a simple polygon optimally, either an optimal parallel algorithm for triangulation 
should be obtained, or a new method that does not need triangulation should be 
developed. 

An optimal parallel algorithm is obtained for triangulating a star shaped polygon 
in Chapter 4. Triangulating any arbitrary polygon in optimal O(logn) time is still 
an open problem. 

The convex hull (of a simple polygon) algorithm proposed in Chapter 5 is optimal. 
It has to be seen whether faster algorithms can be obtained. 

Some fast parallel algorithms for computing visibility polygon are presented in 
Chapter 6. These algorithms work on a powerful CRCW COMMON PRAM. It 
has to be seen whether same bounds can be obtained on a weaker model. Finding a 
0((log log n) 2 ) time algorithm with n processors would be more interesting. If one can 
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avoid the binary search bottle-neck, i.e. avoid binary search that is being performed 
to find intersection of visibility chains, then can be solved in desired resource bounds. 

In Chapter 7 a fully optimal algorithm to get point set from a Voronoi diagram 
is described. It has to be seen whether this algorithm can be extended to the case, 
where points are not in general position. 
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